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ABSTRACT 

Context. The increased sensitivity and high spectral resolution of millimeter telescopes allow the detection of an increasing number of 
isotopically substituted molecules in the interstellar medium. The ratio is difficult to measure directly for carbon containing 

molecules. 

Aims. We want to check the underlying hypothesis that the ^^C/ ^^C ratio of nitriles and isonitriles is equal to the elemental value via 
a chemical time dependent gas phase chemical model. 

Methods. We have built a chemical network containing D, ^^C and molecular species after a careful check of the possible 
fractionation reactions at work in the gas phase. 

Results. Model results obtained for 2 different physical conditions corresponding respectively to a moderately dense cloud in an 
early evolutionary stage and a dense depleted pre-stellar core tend to show that ammonia and its singly deuterated form are somewhat 
enriched in in agreement with observations. The ratio in N 2 H+ is found to be close to the elemental value, in contrast to 

previous models which obtain a significant enrichment, as we found that the fractionation reaction between and has a barrier 
in the entrance channel. The large values of the N 2 HV^^NNH^ and N 2 H^/ ratios derived in L1544 cannot be reproduced in 

our model. Finally we find that nitriles and isonitriles are in fact significantly depleted in questioning previous interpretations of 
observed C^^N, HC^^N and H^^NC abundances from containing isotopologues. 
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1. Introduction 


Understanding isotopic abundances on a large scale is a major 
field of interest which has received a great deal of attention for its 
application to terrestrial environments (ocean, meteorites), the 
solar system (planets, comets) and in galactic interstellar space. 
The variation of isotopic ratios may give us some information 
about the link between solar system objects and g alactic inter¬ 
stellar environments as discussed by |Aleon| ( |2Q10| ). We princi¬ 
pally focus our study on interstellar environments where low 
temperature conditions may significantly impact the isotopic ra¬ 
tios of the molecular content. Isotopic molecules are detected 
in a variety of environments and offer an additional tool to de¬ 
termine physical conditions as they usually do not suffer from 
opacity problems. 


Early modeling studies on and isotopic enrichment 
were performed by [Langer et al.| ( |1984| ) who introduced differ¬ 
ent isotopic exchange reactions, relying on previous theoretical 
and experimental studies by [Watson et al.| 1976} [Smith & Adams 
[1980[ Possible effects of selecti ve photodissociation of CO have 
subsequently been emphasized ([Glassgold et al.|1985t [Le Bour- 


[lot et al.[|l993| [Visser et aL][20b9[ ), which tend to increase the 

CO/^^CO ratio. The use of the C N radical as a tracer of 
isotopic ratio has been raised by [Savage et al.[ ( [2002[ ); [Milam 
et al.[ ( [2005] ), who studied the corresponding gradient as a func¬ 
tion of the galactic distance. The actual value of the iso¬ 


topic ratio in the local Interstellar Medium (ISM) is assumed to 
be 68 ( [Milam et al.|2005[ ). 

The possibility of nitrogen isotopic fractionation in interstel¬ 
lar clouds has been investigated by [Terzieva & Herbst] p000[) 
who suggested various isotopic exchange reactions. Rodgers 


[& Chaml^ ( [20Q4[ [2008a|b[ ) used these suggested reaction rate 


constants to predict nitrogen isotopic fractionation in chemical 
models of dense interstellar molecular cores. They specifically 
discussed the role of the atomic to molecular nitrogen ratio in 
the fractionation process and the possible link between nitro¬ 
gen hydrides and CN containing molecules (nitriles). The cor¬ 
responding observations are sparse however and difficult as the 
elemental ratio is high (441 + 5), as determined by the 

recen t Genesis solar wind sampling measurement ( [Marty et~ar 


[201 1[ ) and assumed to hold in the local ISM. In addition, the zero 
point energy (ZPE) differences involved in nitrogen fractionation 
reactions are small an d the predicted corresp onding chemical en¬ 
richment is moderate. 


Lucas & Liszt 


corresp i 

]( |1998| ) 


reported HC^^N ab¬ 
sorption in diffuse clouds located in front of distant quasars with 
a HCN/HC^^N ratio of 270+ 27, close to value reported for the 
Earth. However, various new observations of isotopic nitrogen 
containing molecules have been reported, including substi- 
tuted amm onia and deuterated ammonia ( Gerin et al.|2009[[Lis[ 
et al.|201()|) , the diazenyl ium ion (N 2 H^) ( Bizzocchi et al.|20f0 


2013 


Daniel et al.|2'0T3l), CN and HCN ([Ikeda et al.|2002[[PiH 


lai et al. |2007t|AdWe & Ziurys|20T2l |Hily-Blant et al.|20i3bl 
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IDaniel et al.|20T3l ), HCN and HNC dWampfler et al.|2014D . The 
strong depletion found in variants of the N 2 H^ isotopologue 
strongly contradicts model predictions ( [Gerin et al.|2009| ), which 
motivates a reinvestigation of the chemical processes at work. 
With this in mind, the link between deuterated chemistry and the 
possible role of ortho/para molecular hydrogen has been studied 
by |Wirstrom et al.| ( |2012| ). 

We analyse in Section the various possible isotopic ex¬ 
change reactions that are involved for carbon and nitrogen con¬ 
taining molecules. Indeed, most nitrogen fractionation observa¬ 
tional results for CN containing molecules involve only and 
species so that the measure of the nitrogen isotopic ratio 
assumes a fixed fraction. We examine and extend the pi¬ 

oneering study of |Terzieva & Herb^ ( |2000| ) and check for the 
possible presence of barriers in the entrance channels of isotopic 
exchange reactions through theoretical calculations. We also up¬ 
date the zero point energy (ZPE) values involved and derive the 
corresponding exothermicity values. We present our new chem¬ 
ical model in Sectionj^and compare with available observations 
and other models. Our conclusions are presented in Section 


2. Chemical reactions involving isotopic substitutes 
of and 

2.1. and exchange reactions 

At very low temperatures, isotopic exchange reactions may only 
occur if no barrier is present between the interacting atoms, ions 
and molecules or if tunnelling plays an important role. Experi¬ 
mental information is crucial and we constrain the evaluation of 
rate constants using that information. If no experimental data are 
available, we apply theoretical methods to determine the pres¬ 
ence of a barrier: A first technique uses DET (Density Eunc- 
tional Theory) calculations (with the hybrid M06-2X functional 
developed by |Zhao & Truhlar| ( |2008| ) which is well suited for 
thermochemical calculations, associated to the cc-pVTZ basis 
set using GAUSSIAN09 software). The alternative is provided 
by the MRCI-rQ method (with the aug-cc-pVTZ basis set). Eor 
barrierless cases, we derive the reaction rate constants by using 
a simple capture theory ( [Georgievskii & Klippenstein[2005| ) for 
both ion-neutral and neutral-neutral reactions. We consider four 
different families of isotopic exchange reactions: 

“ A : direct reactions. The proton transfer in the N 2 H^ -l- 
^ -I- N 2 reaction can serve as an example. 

In this case and for reactions without a barrier, the reac¬ 
tion rate coefficient of the forward reaction is equal to the 
capture rate constant multiplied by a probability factor/(5, 
M), f{B, M) depending on the rotational constant, mass and 
symmetry values of the reactants and products. In reactions 
involving^^N and the mass ratio of reactants and prod¬ 
ucts are very close and M) = CFo^^^^^incQ channel/t^exlt channels 
(The symmetry number a is equal to the number of path¬ 
ways). The reverse reacti on is calculated from the equi¬ 
librium constant K, as in [Terzieva & Herbst| (|2000l): K = 
kfikr = f(B, M) X exp(AE/kT). 

- B : reactions involving adduct formation leading to direct 
products without isomerization. As an example, we refer 
to -h N 2 ^ -h We first assume that the 

high pressure rate constant is equal to the capture rate con¬ 
stant (for reactions without a barrier). We apply statisti¬ 
cal theory for the system at thermal equilibrium so that 
kf kr = kcapture- Erom the equilibrium constant expres- 
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Sion, we then derive kf = kcapmre x [/(g,M)+Lp(-A£/^7-)] “d 

u _ u w exp{-^E|kT) 

- /^capture ^ [f(B,M)+exp(-AE/kT)] 

- C : reactions involving adduct formation with isomerization 

pathways. Such a case holds for + HCN ^ + H^^CN. 

We again assume that the high pressure rate constant is given 
by capture theory (for reactions without a barrier). The iso¬ 
topic isomerization reaction competes with the dissociation 
of the adduct. The rate constant depends on the location of 
the transition state and statistical calculations are generally 
required to estimate the isomerization reaction rate constant. 

- D : other reactive exothermic channels exist. The exchange 
reaction is discarded generally (the possibility of N atom ex¬ 
change in the -r CN and in -r C 2 N reactions is dis¬ 
cussed). 


The knowledge of the exoergicity values AE is also a major con¬ 
cern. They are obtained from the differences of the ZPEs be¬ 
tween products and reactants. We recall in the Appendix the cor¬ 
responding expressions and derive their values by using the most 
recent determinations of spectroscopic constants. 

We summarize in Table the different isotopic exchange re¬ 
actions considered and display the corresponding reaction rate 
constants. Detailed information is provided in the online mate¬ 
rial on the theoretical methods used for the different systems. 
The reactions involving are displayed in the upper part of 
the Table. We also consider isotopic exchange reactions in 
the lower part of Table Table p] shows two main d iscrepancies 

compared to previous studies by |Terzieva & Herb'^ ( |2000| ) : The 
exchange reactions between atomic and N 2 H^, HCNH^ are 
found to be unlikely to occur as significant barriers arise in the 
complex formation step. A similar result is obtained for 
ex change with NO, whereas these reactions had been included 
in [Terzieva & Herbst] ( |2000| ). The exchang e reaction between 


atomic and CN, which was suggested by Rodgers & Charn 
|ley| ( |2008^ , is found to be plausible. Additional possibilities of 
exchange have also been considered such as the reaction between 
^^N and C 2 N. As far as possible fractionation is concerned, 
we find that CN could be enriched in through the exchange 
reactions of CN with ^^C and ^^C^. However, such a mecha¬ 
nism does not hold for HNC a s atomic carbon is found to react 
with HNC ( Loison et al.|2014 ). ^^C enrichment of HCN is also 
found to be unlikely as the calculated transition state lies above 
the entrance level in the hypothetical isomerization process (C 
mechanism). 


2.2. Ammonia synthesis 

Ammonia synthesis proceeds mainly through a chain of reac¬ 
tions starting with ^^N^ and H 2 , as t he reaction between N and 
H 3 ^ has been shown to be inefficient ( Milligan et al.|2000 ). 


2.2.1. The + H 2 reaction and isotopic substitutions 

This almost thermoneutral reaction deserves a special mention 
and has received considerable attention. |Le Bourlot| \\99\) first 
pointed out the possible role of ortho-H 2 in the interstellar 
chemistry of ammonia as the energy of its J=1 rotational level 
almost compensates the small endothermicity of the reaction 


14 


N+ -r H2 


NH+ -h H. 


Dislaire et al. (2012) subsequently re¬ 


analysed the experimental data ( [Marquette et al.|1988] ) and sug 
gested new separate expressions for the re action rate wi t h p-H 2 
and 0 -H 2 . Similar results were obtained by |Zymak et al.| ( [MT3] ), 
who also emphasized the possible role of the fine structure level 
of ^^N^. We follow the prescription derived by Dislaire et al. 
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( | 2012 | ) and extend their analysis to deuterated forms and those 
including , as displayed in TableIn the case of substi¬ 
tuted compounds, we have taken into account the (small) addi¬ 
tional term due to the change in ZPE. These expressions should 
be used with caution as in their kinetic expression, we consider 
that the exponential term represents the enthalpy difference be¬ 
tween the products and reactants. The capture rate constant of 
-h HD reaction is about (2/3)^-^ smaller than the rate of the 
-f H 2 reaction, due to the different mass dependences. The 
formation of is favored at low temperatures. 


2.2.2. NHf + Ho 


The final step of ion-molecule reactions leading to ammonia for¬ 
mation is the reaction between NH 3 and H 2 , giving NH^. Al¬ 
though exothermic, this reaction has a strong temperature depen¬ 
dence, displaying a minimum at T ~ IQOK and a slow increase at 
lower temperatures ( Barlow & Dunn|I987| ), which is interpreted 
by the presence of a barrier to complex formation, as discussed 
in |Herbst et al.| ( |I99T] ). At temperatures close to I OK, the reaction 
is likely to proceed through tunneling which may take place or 
H atom abstraction. However, the NH 3 reaction with D 2 is found 
to be slower when the temperature decreases as tunneling is not 
efficient with deuterium. We thus reconsider the isotopic variants 
of this reaction, as shown in Table where we give the present 
react ion rates compared to prev ious values which were derived 
from |Anicich & Huntress| ( p^ 86 | ). These values are indeed differ¬ 
ent from those used in our previous studies ( [Roueff et al.|2QQ5| ) 
where we assumed the same rate for the channels resulting from 
the reactions of NH 3 and isotopologues with HD, based on pure 
statistical considerations. We discuss the resulting modifications 
in Section O Identical reaction rate coefficients are used for the 
isotopically substituted reactions. 


3. Models 

3.1. General features 


Chemical reactions involving nitrogen atoms and CH, CN and 
OH have been studied experimentally at low temperatures and 
shown to be less efficient than prev iously thought at low temper¬ 
atures ( Daranlot et al. Pol'S)[20 13 j ), which was confirmed by the¬ 
oretical studies. The corresponding reaction rate constants have 
been imple mented in the KIDA chemical data base ( jWakelanij 
|et al.| 2 QI^ and we have updated our chemical network accord¬ 
ingly. We also include the reactions discussed in [Loison et al.j 
( |20I4| ) in their study of HCN / HNC chemistry. We explicitly 
include deuterium, and molecular compounds in our 
chemical network. The reactions displayed in Table have been 
included, which allows us to test the hypothesis of a constant 
isotopic ratio to derive the ratio in con¬ 

taining molecules. We take into account the role of the ortho/para 
ratio of molecular H 2 in the H 2 D^ -l- H 2 and -l- H 2 

chemical reactions in an approximate way: we do not compute 
the full ortho/para equilibrium as in the models of [Flower et al.j 
2006[[Pagani ct al.|20IH[Faure et al.|20I3[ a nd rather introduce 


it as a model parameter, which can be varied. [Faure et al.[ ( [2QI3[ ) 
have found a value of 10“^ for temperatures below 15K. 

Apart from exchange reactions, reactions involving isotopic 
molecules are assumed to have the same total rate constant as 
those involving the main isotope except for the reaction of 
with H 2 /HD/D 2 . The various reaction channels are obtained from 
statistical considerations, in the absence of experimental infor¬ 
mation. We restrict carbon containing molecules to 3 carbon 


atoms, nitrogen containing molecules to 2 nitrogen atoms and 


consider full deuteration as in our previous studies (Roueff et al. 


2005). Within these constraints, the number of species consid¬ 


ered in the model is 307 linked through more than 5400 chemical 
reactions. We consider two different models as displayed in Ta- 


Table 4. Model d efinitions. The C / and N / ratios are re spec- 
tively taken as 68 ( [Milam et al.|2005[ l and 440 ( [Marty et al.|2011| ). 



Model (a) 

Model (b) 

density nn (cm“^) 

2x 10 ^ 

2x 105 

Temperature (K) 

10 

10 

cosmic ionization rate per H 2 (s“^) 

1.3 X 10-1^ 

1.3 X 10-1^ 

He/H 

0.1 

0.1 

C/H 

4.15 X 10-5 

1.4 X 10-5 

N/H 

6.4 X 10-5 

2.1 X 10-5 

0/H 

6 X 10-5 

2.0 X 10-5 

S/H 

8.0 X 10->^ 

8.0 X 10->^ 

Fe/H 

1.5 X 10-^ 

1.5 X 10-^ 


ble|^ Model (a) may be considered as a template of TMCI, and 
assumes a density of hydrogen nuclei vih = 2x cm~^. The el¬ 
emental abundance of carbon relative to hydrogen nuclei is taken 
as 4.15 X 10“^ to reproduce the derived relative abundance of CO 
( [Qhishi et al.|1992[ ). We derive the oxygen elemental abundance 
by imposing a C/0 ratio of 0.7 appropriate for TMCI and take 
the nitrogen elemental abundance used by [Fe Gal et al.[ ( [20I4[ ) 
in their work on nitrogen chemistry. The elemental abundance of 
sulfur is not well constrained and we have taken the low metal 
case value of 8.0 x 10“^. Model (b) is more representative of a 
pre stellar core with a density of 2 x 10^ cm“^ similar to FI34N 
or Barnard 1 (BI) where the elemental abundances of carbon, 
oxygen and nitrogen are reduced by a factor of 3 to account for 
depletion. T = I OK in both cases and the cosmic ionization rate 


^ per H 2 is 1.3 x 10 s ^ as in Le Gal et al. (2014). 


3.2. Results 

We summarize our results obtained with a 10“^ value of the o/p 
ratio of H 2 in Table and give some observational values for 
comparison. Time dependent effects may be visualized from the 
values reported at 10 ^ years and at steady state for model (a). 
Steady state is reached after a few 10^ and 10^ years respectively 
for models (a) and (b). As there are fewer enrichment re¬ 
actions than previously assumed ( Terzieva & Herbst|2000[ ) most 
nitrogen containing species are found to have isotopic abundance 
ratios close t o the solar value = 440) given by the Gen¬ 

esis mission ( [Marty et al.|2QII[ ). 


4. Discussion 

We display the time dependence of various isotopic ratios and 
fractional abundances relative to H 2 and discuss the chemical 
behavior involved in the fractionation processes for the two re¬ 
ported models. We first consider in Figure the reservoirs of 
nitrogen, atomic and molecular nitrogen as well as N 2 H^ ions 
which are chemically linked to N 2 . Atomic nitrogen becomes de¬ 
pleted in after about 10 ^ years whereas molecular nitrogen 
is slightly enriched. These evolution times are also required to 
build significant amounts of molecular compounds which com¬ 
pare satisfactorily with available observations. The overall de¬ 
pendence of N 2 H^ follows closely that of N 2 as it is formed from 
N 2 + Hs^ reaction, with a slight decoupling between ^^NNH^ 
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Fig. 1. Upper panel : time dg)endence of N / isotopic ratios in atomic and molecular nitrogen and ions, (a) and (b) correspond to 

the models defined in Table Q the black heavy dotted line represents the elemental N / Lower panel : ti me dependence of t he fractional 
abundances relative to H 2 of N, N 2 and N 2 H+ for models (a) and (b). The value reported for N 2 H+ towards TMCl ( [Ohishi et al.|1992| ) is displayed 
as a horizontal green dashed line. 


and at long evolution times. We find that the isotopic 

ratio of the N 2 H^ ions displays an almost constant value close to 
the solar value after some 10^ years. They are in good agreement 
with observations in B1 but disagree by a factor of 2 for LI544. 
The trend that is less abundant than is repro¬ 

duced in our results, as a result of the differences of endothermic- 
ity in their reactions with N 2 . We checked that introducing ^^N 2 
and species containing two atoms had no effect on these 
ratios. We could not find any gas-phase mechanism able to gen¬ 
erate such a large ratio in pre-stellar core conditions. The high 
isotopic ratio found in LI544 implies an equivalently large ratio 
for molecular nitrogen, which is in strong contradiction with our 
fin dings. These results are m arkedly different from those derived 
by Hily-Blant et al. ( 2013bt who found a moderate enrich¬ 
ment as these authors introduced the -f N 2 H^ fractionation 
reaction which we have found not to occur. 0 


4.1. Nitrogen hydrides 

We display in Figure]^ the time evolution of the isotopic ratios 
of nitrogen hydrides and of their fractional abundances relative 
to H 2 . NH 3 and NH 2 have a very similar behavior as they both 
result from the reaction chain starting with the -f H 2 re¬ 
actions. The species NH 2 and NH 3 are found to be enriched 
in due to the -i- 0 -H 2 reaction which has a slightly 
smaller (weak) endothermicity, as reported in Table than the 
corresponding -f 0 -H 2 reaction, a difference which slightly 
favors formation. Nevertheless, despite the large uncer¬ 

tainties regarding the rate of the -f H 2 reaction, our results 
are in fair agreement with the results of |Daniel et al.| ( [ 2 QT^ for 


^ These authors also interchanged the endothermicities of the 
and reactions with N and N 2 . 
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prestellar core Bl. As the -r HD ^ ND^ -r H reaction has 
a smaller endothermicity than the -r H 2 ^ NH^ -r H reac¬ 
tion, the modeled ^'^NH 2 D / ^^NH 2 D ratio behaves somewhat 
differently than the ^^NH 3 / ^^NH 3 ratio. The ratio exhibits large 
variations around 10 ^ years and becomes smaller than that of 
NH 3 at large times and at steady state. The observed values are 
compatible with calculations at sufficiently large times. The val¬ 
ues of this ratio reported in |Gerin et al.| ( |2009| ) have been found 
to be too large as a result of assuming a single excitation rota¬ 
tional temperature. The future availability of collisional excita¬ 
tion rates of NH 2 D by H 2 ( [Daniel et al.|2Q14| ) may give rise to 
additional changes. NH does not follow the same trend as NH 2 
and NH 3 and is only slightly enriched in ^^N.As discussed by 
[Hily-Blant et al.|p013b| ), NH is mainly formed from the disso¬ 
ciative recombination (DR) of N 2 H^, a reaction which has been 
recently revisited by [Vigren et al.[ ( [ 20 T^ who derive a branching 
ratio towards NH of 7%. We have checked that this analysis still 
holds for model (b) even if the reaction of N 2 H^ with CO be¬ 
comes more efficient. The formation route of NH through NH^ 
recombination may take over for highly depleted CO. 


4.2. Nitriles and isonitriles 


Deriving ^^N isotopic ratios of CN, HCN and HNC from ob¬ 
servations is a difficult challenge as the transitions of the main 
isotopologues are optically thick. Then, most of the reported ob¬ 
servational values of the ^"^N / ^^N molecular ratios are obtained 
from the ratios of the minor isotopologues ^^CN / C^^N, H^^CN 
/ HC^^N and HN^^C / H^^NC which is subsequent ly multiplied 
by an assumed C / ^^C value, usually taken as 68 ([Milam et al. 
[2^5] ). 
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Fig. 2. Upper panel : Time dependence of N / isotopic ratios in nitrogen hydrides, (a) and (b) correspond to the models defined in Table 
The black heavy dotted line represents the elemental N / Lower panel : time dependence o f the fractio nal abundances re lative to H 2 of NH, 
NH 2 , NH 3 and NH 2 D for models (a) and (b). The values reported for NH 3 ( [Ohishi et al.|1992| ) and NH 2 D ( [Tine et al.|2000| ) towards TMCl are 
displayed as horizontal blue and red dashed lines respectively in the left panel. 




Fig. 3. Time evolution of ratios in CN, HCN and HNC for models (a) and (b). The elemental ratio is displayed as heavy dotted 

line. 


4.2.1. ratios 


We test these hypotheses in our models by explicitly introducing 
fractionation reactions of as discussed in Section]^ Table 
shows that the isotopic ratios of CN, HCN and HNC vary 
both with time and density. Figuredisplays the ^^CN / C^^N, 
H^^CN / HC^^N and HN^^C / H^^NC ratios as a function of time 
for the two considered models. The deviation from the elemental 
ratio of 6.48 is significant for HCN, HNC and CN. 


4.2.2. chemistry 

We now consider the time dependences of the ^^C / ^^C isotopic 
ratios in CN, HCN and HNC species as displayed in Figure 
The time dependent ratio displays large variations, which is due 
to the fact that there are various reactions incorporating ^^C in 
the molecules. The elemental value of the ratio (68) is fulfilled 
in a narrow range around 1 Myr for model (a) and 2 x 10^ yrs 
for model (b) but steady state values are significantly different 
except for CN. This relatively complex behavior results from the 
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Fig. 4. time dependence of C / isotopic ratios in 


CN, HCN and HNC. (a) and (b) correspond to the models defined in Table|^ 






Fig. 5. Upper panel: Time dependence of C / isotopic ratios in C, CH, CO and HCO^. The black heavy dotted line represents the elemental 
Lower panel: Time dependence of the fractional abundances relative to H 2 of C, CO, CH and HCO^, for models (a) and (b). The 
observational values towards TMCl ( [Ohishi et aL|1992| ) are displayed as horizontal dashed lines with corresponding colors in the left panel, (a) 
and (b) correspond to the models defined in Table 


many different reaction channels involved in chemistry. We 
then also consider other containing species and display the 
ratio in C, CH, CO and HCO^ in Figure as well as 
their fractional abundances relative to H 2 . The ^^C/^C isotopic 
ratios of the various molecules are highly dependent on the evo¬ 
lution time. The transition from gas-phase atomic carbon toward 
CO controls the ^^C enrichment. As long as there is still a rel¬ 
atively high carbon atom concentration in the gas phase, there 
is enough free ^^C to allow strong enrichment of CN through 
the ^^C -f CN reaction. When CO molecules become the reser¬ 
voir of carbon, even if in that case the ^^C concentration is low. 

Article number, page 6 of[^ 


the ^^C^ -f ^^CO ^^C^ -f ^^CO reaction still leads to a smal l 
^^CO enrichment ( Langer & Penzias|l990 Milam et ^|2005 ) 


Although this small excess is not measurable in CO, significant 
amounts of ^^C are locked up in CO and most of the other carbon 
containing species become depleted in ^^C, as found for CH and 
other carbon chains. This effect is seen in Figure for the iso¬ 
topic ratio of CH in model (a) at steady state where CO is slightly 
enriched, leading to a significant depletion of ^^C in CH. Our re- 
sults for C and HCO^ are similar to those of Fumy a et al. ( [2011 1 
who studied the ^^C fractionation of multiple carbon chains by 
explicitly introducing the dependence of the ^^C position in the 
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chain. We see that HCO^ is marginally enriched in at steady 
state whereas HCN and HNC are significantly depleted in 
We also note that CN and HCO^ may react via proton transfer 
to give CO + HNC^ and have checked , at the DFT level, the 
absence of any barrier. The corresponding reaction rate is 2.2 x 
10“^ cm^ s“^ when using the capture rate theory. This 

reaction has not been included in any chemical database up to 
now. 


HNC has been shown recently to react with atomic carbon 
( Loison et al.|2014| ), which leads to the different steady state iso¬ 
topic ratios obtained for HCN and HNC. The CN chemistry is 
then somewhat decoupled from that of HCN and HNC. HCN 
and HNC are formed at relatively long times via CN -f Hs^ ^ 
HNC+ /HCN+ -f H 2 , followed immediately by H NC+ -r H 2 ^ 
HCNH^ -I- H giving back HCN and HNC via DR ( [Mendes et al. 


2012). With the adopted elemental abundances, the main CN de¬ 


struction reactions are however O -h CN and N -f CN so that the 
HCNH+/HCN+/HNC+/HCN/HNC/CN network is not a closed 
system in contrast to models including coupled gas-grain chem¬ 
istry ( [Loison et al.|2014| ). 


4.2.3. fractionation 

We now display also the N / ^^N fractionation in nitriles as well 
as that of NO in Figure The time dependences of the isotopic 
ratios are markedly different in the two models, except for NO. 
Whereas CN, HCN and HNC are somewhat enriched in ^^N for 
model (a) both at intermediates times and steady state, the oppo¬ 
site result is obtained in model (b) after 10^ years. We display as 
well the time dependences of the fractional abundances of these 
molecules in Figure in order to better understand the previ¬ 
ously described differences observed in the various fractionation 
ratios. 

Our model (a) is intended to be representative of TMCl, a 
moderately dense cloud in an early evolutionary stage. We see 
that the various abundances and ratios are very sensitive to the 
chosen "age", assumed to be the relevant chemical evolution 
time. Considering a TMCl age of IMyr leads to a reasonable 
agreement with the sparse observations available, see Table |5 
Model (b) is more likely to be representative of a denser evolvec 
molecular cloud such as L134N, L1544, Barnard Bl, ... where 
elemental C, O, N are partially depleted through sticking on 
grains. Steady state values obtained with significant depletion 
conditions ( [Roueff et al.||2005] ) may be used as corresponding 
proxies. We see that ammonia isotopologues are satisfactorily 
accounted for in our model. The agreement between observed 
and calculated ^^C/ ^^N ratios is somewhat misleading as the 
modeled value mainly results from a significant depletion in the 
^^C species. The ratios involving N / ^^N are found close to the 
elemental value in our model (b) at steady state (although NH 2 
and NH 3 are somewhat enriched in ^^N through the ^^N^ -f H 2 
reaction as explained in subsection [4.1[ ), which is in agreement 
with the fact that no significant gas phase fractionation mecha¬ 
nisms have been found. The occurrence of small ratios in B1 ob¬ 
servations ( [Daniel et al.|2Q13[ ), if real, implies other mechanisms 
at work. An obvious suggestion lies in the processes involved in 
adsorption/desorption on grains and possible surface reactions. 


4.3. Role of the H 2 ortho to para ratio 

We now test the role of the value of the ortho/para ratio of 
molecular hydrogen (OPR) in the context of ammonia chemistry 
and fractionation determination. We run two additional models 


in the frame of models (a) and (b) by changing the OPR by a 
factor of 10 both upwards and downwards. The abundances of 
P-H 2 and 0 -H 2 are expressed respectively as n(p-H 2 ) = 
n(H 2 ) and n(o-H 2 ) = n(H 2 ). This ratio is introduced, 

in addition to the reaction of H 2 with ^"^N^ (^^N^), in the 
reverse reaction of the fractionation of H 3 ^, namely the H 2 D^-l- 
H 2 reaction which plays a signi ficant role in the deuterium 
fractionation of various molecules ( [Pagani et al.|2Qll[ ) as shown 
below: 

H2D++P-H2 ^ H3 +HD ,ki(cm 3 s-i) 
H2D++0-H2 ^ H3++HD ,k2(cm3s-i) 

with ki = 2.0 X 10“^ X exp(-232/T) and k 2 = 2.0 x 10“^ x 
exp(-61.5/T). 

Table 7. Dependence of the fractionation ratios on the o/p ratio of H 2 
for model (b) at steady state. 



OPR = 10“^ 

OPR = 10“^ 

OPR = 10“^ 

electronic fraction 

2 . 2 x 10 -'' 

1.7 X 10“^ 

1.1 X lO-"" 

N/15n 

457 

455 

445 

2 X N 2 / ‘^NN 

436 

437 

438 

NH / 

425 

421 

418 

NH/ND 

9 

9 

11 

NH 3 / 

396 

387 

416 

NH2D/ 15NH2D 

322 

336 

376 

NH 3 / NH 2 D 

II 

18 

33 

N2HV N‘ 5 NH+ 

425 

423 

417 

N 2 H-"/ 15 NNH+ 

434 

433 

433 

N2HV N2D+ 

8.7 

8.6 

9.7 

CN / '^CN 

65 

63 

60 

CN/C'^N 

449 

445 

438 

HCN / H'^CN 

II5 

114 

107 

HCN/HC'^N 

467 

453 

458 

HCN / DCN 

25 

22 

18 

HNC/ HN^^C 

120 

121 

118 

HNC/ H^^NC 

461 

446 

453 

HNC/ DNC 

19 

16 

12 

HCO+ / H‘3C0+ 

58 

56 

51 

HCO+/DCO+ 

8.6 

8.4 

9.4 


Tables and [7] display the fractionation values for the three 
different assumed values of the OPR for models (a) and (b) and 
Figures and display the corresponding time evolutions. We 
see that the curves are almost superposable for times less than 1 
Myr for model (a) and less than several 10^ yrs for model (b). 
The variations are significative at steady state. We find that the 
competition between the destruction channels of ^^N^ through 
its reactions with 0 -H 2 and CO plays a major role. If 0 -H 2 is 
the most efficient destruction channel, which occurs typically for 
OPR > 200 X xco at I OK, where xco is the fractional abundance 
of CO relative to H 2 , formation of NH^ proceeds efficiently. In 
the opposite case ^^N^ is mainly destroyed through reaction with 
CO to yield CO^, leading rapidly to HCO^, and NO^ (which 
does not react with H 2 ). As the DR rate coefficients of poly¬ 
atomic ions increase significantly for larger polyatomic ions (the 
NH 4 -r e“ reaction is 3 times more rapid than HCO^ -r e“ and 
20 times more rapid than H 3 ^ -i- e“), these two different chan¬ 
nels impact the electron abundances and then affect many other 
species. Moreover, the reaction between NH3 and H3^ leads to 
NH4 -I- H 2 . NH4 then reacts mainly with electrons, recycling 
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Fig. 6. Time dependence of N/ isotopic ratios in 


CN, HCN, HNC and NO. (a) and (b) correspond to the models defined in Table|^ 


back to NH 3 which subsequently converts to NH|, amplifying 
the electron loss through DR. These effects occur when CO be¬ 
comes the main carbon reservoir, for sufficiently large evolution 
times. As a result, the electron fraction is found to decrease when 
the OPR increases at large evolution times and at steady state. 
The decreasing electron abundance with increasing OPR values 
acts to diminish the effect of all DR reactions and as a result to 
increase the abundances of and H 2 D^ since DR is the main 
destruction channel in our conditions. The abundance of H 2 D^ 
is even more enhanced through the -l- HD reaction. Then the 
abundances of other deuterated molecular ions produced through 
deuteron transfer reactions of abundant neutrals with H 2 D^ are 
increased as well, which impacts the subsequent formation of 
neutral deuterated compounds produced in the DR reactions. 
Eventually, the HCN/DCN and HNC/DNC ratios are shown to 
decrease for increasing OPR values. This example demonstrates 
the complexity of the interplay between the different chemistries. 

5. Conclusions 

We have built for the first time an isotopically substituted gas 
phase chemical network including D, and species and 
included it in a time dependent chemical model. Our model is 
based on a careful analysis of the possible gas phase mecha¬ 
nisms involved in carbon and nitrogen fractionation by scrutiniz¬ 
ing the few available experimental studies and performing DFT 
and ab-initio quantum calculations on hypothetical reactions to 
check the possible presence of barriers in the reaction channels. 
One important result obtained is that the fractionation reaction 
of with N 2 H+ is unlikely, due to the pres ence of a barrier, in 
contrast to the previous hypothesis made by [Terzieva & Herbst 
P000| ). As a result, the modeled isotopic ratios involved in the 
isotopologues of N 2 H^are found to be very close to the elemen¬ 
tal values and are similar to each other, in contradiction with 
observations towards LI544 ( [Bizzocchi et al.|2()T3] ). The avail¬ 
ability of new collisional rate coefficients for the N 2 H^- H 2 sys- 
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tern ( [Lique et al.|2015| ) may however modify these conclusions. 
We also discarded the -r HCNH^ and -r NO exchange 
reactions through similar arguments. Tentative reaction rate co¬ 
efficients are also proposed for carbon fractionation reactions in¬ 
volving and with CN and C 2 . Additionally, we have 
explicitly considered the various isotopologues involved in N^ 
-r H 2 reaction, assuming that the energy defect involved in the 
reactions of ^"^N^ with para-H 2 is a "real" endothermicity. This 
leads to a slight decline of the exponential term when ^^N^ re¬ 
acts with H 2 and with HD compared to ^"^N^. This explains sat¬ 
isfactorily that ^^NH 2 D is found to be more enriched in ^^N than 
^^NH 3 in the observations. Comparison between observations of 
nitriles and isonitriles and simulated values is much more ques¬ 
tionable, as carbon and nitrogen chemistries are interdependent. 
Observed isotopic ratios are usually large and suffer from large 
error bars due to opacity effects in the main isotopologue and 
difficulties linked to nuclear spin effects. Whereas the various 
isotopologues follow a similar evolution, the isotopic ratios dis¬ 
play significant variations due to slight shifts in the position of 
the maximum fractional abundances. A reasonable agreement is 
obtained between the observed ^^C / ^^N ratios for most of the 
species in L134N and Barnard B1 and steady state model val¬ 
ues. Our model results show a strong depletion in ^^C and a 
near elemental ratio for ^"^N / ^^N, whereas observations are usu¬ 
ally interpreted by assuming an elemental ratio for ^^C contain¬ 
ing species which leads to the incorrect assumption of ^^N en¬ 
richment. These considerations are undeniably dependent on the 
chosen elemental abundances, and in particular to the assumed 
C/0 ratio. We additionally point out a somewhat unexpected ef¬ 
fect of the ortho to para ratio of H 2 , which affects significantly 
the fractional ionization, and consequently the level of deuterium 
fractionation through the respective orders of magnitude of DR 
rate coefficients of polyatomic molecular ions. The importance 
of coupling C, O and N chemistries is emphasized. 
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Fig. 7. Time dependence of CN, HCN, HNC and their isotopologues. (a) and (b) correspond to the models defined in Table]^ 
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Table 1. Isotopic exchange reactions. 


label / 
comment 


Reaction 

k/‘ 

(cm^ s 


AE* 

(K) 

(1)A 

N“N + N 2 H+ 


N“NH+ + N 2 

2.3 X 10-“ 

0.5 

10.3 

(2) A 

N“N + N 2 H+ 


“NNH+ + N 2 

2.3 X 10-“ 

0.5 

2.1 

(3) A 

N“N + “NNff 

h 

N“NH+ +N“N 

4.6 X 10-“ 

1 

8.1 

(4)B 

“N+ + N 2 


14n+ +n“N 

4.8 xl0-“ X =-— 

2+exp{-2?>3IT) 

2 

28.3 

(5)C 

“N + CNC+ 


C“NC+ + “N 

3.8xlO-“x(^)-i 

1 

38.1 

(6)D 

15n+ + 14^0 


14n+ + 15^0 

no react 

- 

24.3 

(7) barrier 

“N + N 2 H+ 


“N + N“NH+ 

no react 

- 

38.5 

(8) barrier 

“N + N 2 H+ 


14n + 15NNH+ 

no react 

- 

30.4 

(9) barrier 

“NNH+ + H 


H + N“NH+ 

no react 

- 

8.1 

(10) barrier 

“N + HCNH+ 


“N + HC“NH+ 

no react 

- 

37.1 

(11)D 

“N + CN 


“N+C“N 

upper limit: 2.0 x 10“^^ x 

1 

22.9 








(12) B 

“N + C 2 N 


“N + C2“N 

1.6xl0-“x(7’/300)“®x 

1 

26.7 





i 

l+ejcp(-26.7/T) 



(13) D 

“N + ‘^NO 


“N + “NO 

- 

- 

24.3 

(14) B 

+ CO 


“C+ + “CO 

6.6 xl0-“ X (7’/300)-'J'45 

1 

34.7 





X exp(-6.5/T) X u) 



(15) A 

“CO + HCO+ 


CO + H“CO+ 

2.6 xl0-“ X (7’/300)-‘''4 

1 

17.4 

(16) B 

“C+ + CN 


“C+ + “CN 

3.82x10-® x(r/300)-0'4 

1 

31.1 





V ^ 

^ l+ejcp(-31.1/r) 



(17) B 

“C + CN 


“C+“CN 

3.0 XlO X l/J) 

1 

31.1 

(18) C 

“C + HCN 


“C + H“CN 

no react 

- 

48.4 

(19) B 

“C + C 2 


“C+“CC 

3.0 XlO X 2+expi-26AIT) 

2 

26.4 

(19) barrier 

“CH + CO 


“CO + CH 

no react 

- 

28.6 


* k/ is the forward reaction rate coefficient. The reverse reaction rate coefficient, , is obtained by k^ = exp(- AE/T). 
Table 2. Reaction rate coefficients of + H 2 and isotopic variants. 


Reaction 


k (cm^ s ') 


Comment 


H 2 

^ NH+ + H 

8.35 X 10-“ X exp(-168.5/T) 

Dislaire et al. 

2012 

H 2 

^ NH+ + H 

4.2 X 10-“ X (T/300)-o “ x exp(-44.5/T) 

Dislaire et al. 

2012 

H 2 

^ “NH+ + H 

8.35 X 10-“x exp(-164.3/T) 

see text 


H 2 

^ “NH+ + H 

4.2 X 10-“ X (T/300)-o “ x exp(-39.7/T) 

see text 



+ 

14 ^+ + 
15n+ + 
15n+ + 

+ HD 
+ HD 
+ HD 
+ HD 

+ D2 
+ D, 


ND+ +H 
NH+ +D 

*5nD+ + H 
*5nH+ + D 
ND+ +D 
‘5nD+ + D 


3.17 X 10-“ X exp(-16.3/T) 
3.17 X 10-“ X exp(-594.3/T) 
3.17 X 10-“ X exp(-9.3/T) 
3.17 X 10-“ X exp(-589.5/T) 
2.37 X 10-“ X exp(-197.9/T) 
2.37 X 10-“ X exp(-190.9/T) 


Marquette et al. 1988 


see text 
see text 
see text 


Marquette et al. 

1988 

Marquette et al. 

1988 
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Table 3. Reaction rate coefficients of NHf + H 2 and isotopic variants at T = lOK. 


Reaction 

k (cm^ 

present work 

s-i) 

old value (*) 

NH+ + H 2 

^ NH+ + H 

8.2 X 10-13 

2.4 X 10-13 

NH+ + HD 

^ NH| + D 

8.2 X 10-13 

1.2 X 10-13 

NH+ + HD 

^ NH 3 D+ + H 

1.0 X 10-13 

1.2 X 10-13 

NH+ + D 2 

^ NH 3 D+ + D 

1.0 X 10-13 

2.4 X 10-13 

NH 2 D+ + H 2 

^ NH 3 D+ + H 

8.2 X 10-13 

2.4 X 10-13 

NH 2 D+ + HD 

^ NH 3 D+ + D 

8.2 X 10-13 

1.2 X 10-13 

NH 2 D+ + HD 

^ NH 2 D+ + H 

1.0 X 10-13 

1.2 X 10-13 

NH 2 D+ + D 2 

^ NH D+ + H 

1.0 X 10-13 

2.4 X 10-13 


NHD+ + H 2 
NHD+ + HD 
NHD+ + HD 
NHD+ + D 2 
ND+ + H 2 
ND+ + HD 
ND+ + HD 
ND^ + D2 


NH 2 D+ + H 
NH 2 D+ + D 

NH D+ + H 
NH D+ + D 
NH D+ + H 
NH D+ + D 
ND| +H 
ND+ +D 


8.2 X 10-‘3 

8.2 X 10-‘3 
1.0 X 10-‘3 
1.0 X 10-*3 

8.2 X 10-13 

8.2 X 10-13 
1.0 X 10-13 
1.0 X 10-13 


2.4 X 10-13 

1.2 X 10-13 

1.2 X 10-13 

2.4 X 10-13 

2.4 X 10-13 

1.2 X 10-13 

1.2 X 10-13 

2.4 X 10-13 


(*) |RouefF et al.| ( [2005| i 
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Table 5. Model (a) and (b) results and observations, ss means stationary state. The value of the o/p ratio of H 2 is taken as 10“^. 



Model (a) 
t= 10^ yrs ss 

Model (b) 

ss 

TMCl 

LI 544 

B1 

electronic fraction 

1.4 X 10-'' 

2.8 X 10-" 

1.7 X 10-" 




N/ 15 n 

440 

456 

455 




2 X N2 / '^NN 

438 

431 

437 




NH / ‘^NH 

429 

428 

421 




NH/ND 

16 

31 

9 




NH3/H2 

6.7 10“® 

1 . 310 -® 

6.0 10-9 

2 10 “^ 



NH3/ 

333 

386 

387 



300*55 (1) 

NH2D/ H2 

3.8 lO-'o 

5.8 10-" 

3.3 10-10 

4 10-10 (9) 



NH2D/ '5NH2D 

215 

276 

336 



230+105 (1) 

NH3/ NH2D 

18 

22 

18 

50 



N2H+/ H2 

4.8 lO-'o 

1.3 10 -‘o 

2.1 10-10 

5 10 - 10 (8) 



N2HV N‘ 5 nH+ 

431 

430 

423 


1050±220 (2) 

400+100 ( 1 ) 

N2H+/ ‘ 5 nNH+ 

437 

432 

433 


1110^240(2) 

> 60019 

N2H+/ N2D+ 

16 

29 

8.6 

12.5 W 


2.9 19 

CN/H2 

6.8 10 “® 

5.5 10-9 

1.2 10-9 

0 

1 

00 



CN/'^CN 

67 

84 

63 



50*[® 19 

CN/C'^N 

430 

449 

445 



240!^f (9 

‘^CN/C'^N 

6.4 

5.3 

7.0 


7.5*1 ( 3 ) 


HCN / H2 

7.4 10-® 

5.9 10-10 

5.4 10-10 

00 

1 

0 

(N 



HCN / H'^CN 

93 

168 

114 



30*2 ( 1 ) 

HCN / HC'^N 

398 

445 

453 



165!50 ( 1 ) 

H'^CN/ HC'^N 

4.3 

2.6 

4.0 

2 - 4.5 (5) 


5.5 +1 19 

HCN / DCN 

43 

96 

22 



20 *®o 19 

HNC/H2 

5.6 10-® 

7.4 10-10 

8.4 10-10 

2 10 "* 



HNC/ HN'^C 

93 

180 

121 

54 - 72 W 


20*5 ( 1 ) 

HNC/ H'^NC 

405 

442 

446 

250 - 330 (4) 


75+25 (1) 
'^-15 

HN'^C / H'^NC 

2.5 

1.75 

3.7 

4.6 + 0.6 W 


3.7 + 1 19 

HNC/ DNC 

23 

66 

16 



2 9+1.1 ( 1 ) 

NO/H2 

1 . 10 “^ 

3.1 10-^ 

4.1 10-^ 

2.7 10“^ 



NO / ‘^NO 

438 

451 

446 




CO/H2 

8.1 10-5 

8.0 10-5 

2.8 10-5 

8 10-5 



CO/'^CO 

68 

67.4 

68 




CH/H2 

1.3 10“® 

1.7 10-9 

1.7 10-10 

2 10 “^ 



CH/'^CH 

74 

154 

71 

>71 



HCO+/H2 

2.5 10-® 

3.610-10 

5.7 10-10 

8 10 -® 



HCO+ / H' 3 C 0 + 

56 

65 

56 



59 (9 

HCO+/DCO+ 

15 

29 

8.4 

50 W 




Referen ces. (l)IDaniel et al.|([M^; (2) [Bizzocchi et a^ ( [20T^ ; (3) |Hily-BlanteraLl(|2013b|); (4)|Liszt & ZiurysldM^ ; (5) IHily-Blant et al 
(2013a|; (6)|Sakai et al.|(|2013|), (7)|Hirano & Liu|(|2014|) assuming CO/ C^®0=500, (8)fohishi et al.|{ 1992|), (9)|Tine et al.|(|2000|), (10)|Gerin et al 
( T9^ 
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Table 6. Dependence of the fractionation ratios on the o/p ratio of H 2 for model (a), ss means stationary state. 



OPR 
t= 10^ yrs 

= 10“^ 

ss 

OPR 
t= 10^ yrs 

= 10-^ 

ss 

OPR 
t= 10^ yrs 

= 10-^ 

ss 

electronic fraction 

5.1 X 10-*^ 

2.2 X 10-' 

4.8 X 10“^ 

2.1 X 10-/ 

3.7 X lO-"" 

4.2 X 10“^ 

N/'^N 

440 

456 

440 

456 

440 

452 

2 X N 2 / ‘^NN 

438 

430 

438 

431 

438 

437 

NH / '^NH 

431 

429 

429 

428 

426 

418 

NH/ND 

16 

31 

16 

31 

18 

23 

NH 3 / 

345 

409 

333 

386 

374 

395 

NH 2 D/ ' 5 NH 2 D 

206 

267 

215 

276 

265 

303 

NH 3 / NH 2 D 

8 

14 

17 

22 

43 

63 

N 2 H+/ N'^NH^ 

431 

429 

431 

430 

429 

421 

N 2 H+/ '5nnh+ 

437 

432 

437 

432 

436 

433 

N 2 H+/ N 2 D+ 

16 

30 

16 

29 

17 

21 

CN / ‘^CN 

67 

85 

67 

84 

67 

70 

CN/C'^N 

432 

449 

430 

449 

429 

434 

HCN / H'^CN 

92 

168 

93 

168 

93 

165 

HCN / HC'^N 

401 

453 

398 

445 

400 

413 

HCN/DCN 

44 

95 

43 

96 

40 

55 

HNC/ HN'^C 

91 

178 

93 

180 

97 

195 

HNC/ H'^NC 

410 

451 

405 

442 

405 

410 

HNC/ DNC 

23 

66 

23 

66 

22 

32 

HCO-" / H'^CQ-" 

57 

66 

56 

65 

54 

55 

HCO-" / DCO-" 

15 

30 

15 

29 

16 

19 
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Appendix A: Exchange reactions 

We discuss below the reactions displayed in Table 


+ N2H+ ^ + N2 

+ N2H+ ^ + N2 

These reactions have been studied experimentally by [Adams 


& Smithj (|1981|) at 80K without differentiating between 
N^^Nlr and^NH+. The total rate is 4.6 x cm^ s”!. 
We thus take half this value for the forward reaction rate con¬ 
stant considering the symmetry factor and we assume no bar¬ 
rier for these rapid reactions. 

We also introduce the ^ -r 

reaction for completeness. 

+ N 2 ^ 

This reaction involves adduct formation. The numerical con- 
stant is computed to reproduce the experimental value of 
Anicich et al.| ( |1977] ) at room temperature. 


- + CNC+ ^ 


T5i^ 


No information is available for this reaction.The bimolec- 
ular exit channels -r NCN and C 2 ^/ N 2 are both en¬ 
dothermic by 170 and 76 kJ/mol respectively. We performed 
DFT calculations (at the M06-2X/cc-pVTZ level) to explore 
the possibility of isotopic exchange. Direct N exchange is 
impossible as it would require simultaneous bond forma¬ 
tion, rearrangement and bond breaking. However, isotopic 
exchange could take place through adduct formation. No bar¬ 
rier was found in the entrance channel for NCNC^ formation 
(exothermic by 167 kJ/mol). The most favorable NCNC^ 
isomerization pathway is cyclic. However the position of the 
c-NC(N)C^ transition state (TS) is highly uncertain and is 
found to be below the entrance channel at M06-2X/cc-pVQZ 
level (-14 kJ/mol) but slightly above the entrance level at 
the RCCSD(T)-F12/aug-cc-pVTZ level (-hl9 kJ/mol). In any 
case, the TS position is close to the entrance level so that the 
NCNC^ back dissociation is favored at room temperature. 
Isotopic exchange may however be enhanced at low temper¬ 
ature and may become the main exit channel. Considering a 
T“^ temperature dependence of the adduct lifetime we tenta¬ 
tively suggest k/ = 3.8 x 10“^^ x (7/300)“^ cm^s“^ 

- 15n+ + 14]V^Q ^ 14]V+ + 15 ^0 
[Terzieva & Herbst] ( |2000| ) estimated the value of this rate 
coefficient from the difference between the Langevin rate 
and the other measured exothermic reactions reported below. 


15 


N+-r 


i^N-r NO+ (ki) 
i^NN^-r O (k2) 
N^^N-r 0+ 
14 n+. 


lS.Hr =-509 kJ/mol 
IsHr =-213 kJ/mol 
^Hr =-402 kJ/mol 
"-r i^NO ^Hr =-0.2kJ/mol 

ki = 5.0 X 10“^^ cm^ s"^ and k 2 = 5.0 x 10“^^ cm^ s“^. 
No data are available for the other two reactions. We find 
a barrier for NON^ adduct formation (261kJ/mol at the 
M06-2X/cc-pVTZ level) and neglect this exchange reaction. 

-r N 2 H+ ^ + n15NH+/15NNH+ 

No information is available for these reactions, which were 
reported as crucial by [Rodgers & Chamley| ( |2008b) ). We per¬ 
formed DFT calculations (at the M06-2X/cc-pVTZ level) 
which showed that isotopic exchange through the addi¬ 
tion elimination mechanism is impossible as the NN(H)N^ 
and NNNH^ ions are metastable, respectivel y 73 and 332 
kJ/mole above the N 2 H^energy. Figure A.l displays 

the various possibilities. In addition, direct nitrogen ex¬ 
change through NN(H)...N^ ^ N...N(H)N^ presents a bar¬ 
rier equal to -r87 kJ/mol. We did not find any NsH^ config¬ 


uration with an energy lower than that of the reactants. We 
then neglect this reaction. 

-r H ^ H -r 

Following a similar suggestion for ex plaining abundanc e 


anomalies of the species of CCH (Sakai et al. 20101, 


we consider the above reaction for which no information 
is available. We performed various DFT calculations at the 
M06-2X/cc-pVTZ level. The HNNH+ ion is 109 kJ/mole 
more stable than H -r N 2 H^ but the addition reaction shows 
a barrier of -l-12kJ/mole. The intramolecular isomerization 
^^NNH^^ N^^NH^ displays also a barrier of -r 172 kJ/mole. 
Then direct isomerization through tunneling is expected to 
be very slow and we neglect that reaction R 
-h HCNH+ ^ -h HC^^NH^ 


This reaction was also suggested by [Terzieva & Herbst 
( [2000[ ) despite the lack of any relevant information. We per¬ 
formed DFT calculations to explore the possibility of iso¬ 
topic exchange as shown in Figure [A.2[ Direct N exchange 
is found to be impossible as it requires simultaneous bond 
formation, rearrangement and bond breaking so that isotopic 
exchange involves adduct formation. Attack by atomic nitro¬ 
gen on either side of the HCNH^ molecular ion leads to a 
metastable system through high energy transition states. We 
did not consider transition states leading to N addition on the 
C=N bond of HCNH^ nor did we examine N insertion into 
the N-H or C-H bonds, as all these pathways lead to species 
located above the reagent energy level. We did not find any 
CN 2 H 2 ^ species lying below the reagent energy level and we 
consider that this reaction cannot take place. 

- -r CN ^ -h 


The ^C^Su) CN (^S^) reaction leads to surfaces in 
Coov symmetry and A" surfaces in Cs symmetry. The quin¬ 
tuplet surface is repulsive at the MRCI-i-Q/aug-ccpVTZ level 
and the ^NCN intermediate is above the N - 1 - CN level. Con¬ 
sidering only the triplet surface, the only barrierless reac¬ 
tion is attack on the carbon atom leading to the groun d state 

|Ma et al.|20T2li 


^NCN intermediate (Daranlot et al. 


2012 


The main exit channel is C -r N 2 after isomerization of the 
NCN intermediate through a tight TS and then back dissoci¬ 
ation may be important and isotope exchange possible. In the 
nominal model we neglect this reaction but some tests have 
been performed to estimate its potential role. The upper limit 
of the isotope exchange rate constant is equal to the capture 
rate constant minus the N -r CN ^ C -r N 2 rate constant. 
That value is notably smaller than the capture one at low 
temperature. We thus propose the upper limit value of 2.0 x 
10-‘o X (T/300)i^® ' 


l+expi-22.9IT) 


constant. 

- C2N^ C2^^N 

The NC^S) -r C 2 N (^H) reaction leads to ^’^H surfaces in Coov 
symmetry and A' -r A" surfaces in Cs symmetry. Figure 
[A.3[ displays the positions of the state energies linked to N 
-r C 2 N which are calculated at the MRCI-hQ/aug-cc-pVTZ 
level (10 e- in 10 OM with the geometry fully optimized at 
the CASCCF level or (14 e- in 14 OM for the CASCCF and 
14e- in 12 OM for the MRCI calculation with non relaxed ge¬ 
ometry). The triplet surface pathway leads to a NCCN adduct 
in a triplet state, corresponding to an excited state of the very 
stable NCCN linear molecule, located -393 kJ/mol below the 
NC^S) -r C 2 N (^H) level. No barrier is present in the entrance 
valley so that this NCCN adduct very likely leads to CN (in 


^ we note in addition that N 2 H^ formation from N 2 + H 3 ^ is exother¬ 
mic by only 72 kJ/mole, so that any produced will not isomer- 

ize into and vice versa. 
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Fig. A.l. Energy diagram of the NNNH+ system. 
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Fig. A.2. Energy diagram of the HCNNH^ system. 


a doublet state) + CN (in a doublet state) products. The oc¬ 
currence of a small exit barrier cannot be excluded but its 
energy should be much lower than that of the ^'^N-l-C 2 en¬ 
trance channel. Triplet surfaces thus cannot lead to isotopic 
exchanges. However the quintuplet surface deserves specific 


attention as the NCCN(^S^) adduct is found at an energy of 
-183 kJ/mol below that of NC^S) -f C 2 N (^H) level (at the 
MRCI-fQ, RCCSD(T) and DFT level of calculations). How¬ 
ever, no exothermic bimolecular exit channel is available on 
this quintuplet surface. We thus conclude that the quintuplet 
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expels*) + CN(*) 


_CN(X^Z") + CN(A^n) 


CNfX^I'^ + CN(X^X'^ 


c, + 'n. 


CCXN 


Fig. A.3. Energy diagram of the CCNN system. 


surface could lead to isotopic exchange and consider finally 
the two following possibilities for this reaction involving the 
very reactive C 2 N radical ( |Wang et ar]|2005[ |2006| ), i.e. the 
triplet channel: 

+ C2N(^n) ^ *5 nccn(^2;'",^ ir) (a.i) 

^ c ‘5 n(22 ;+,2 n) + CN(^ 2 ;+n) (a. 2 ) 


and the quintuplet channel: 

‘5N(''5) + C2N(2n) ^ '5NCCN(^2;+) 

15> 


15x7/4 


NC5) + C2Nrn) 


NC5)+C2Nrn) ^ ^^NCCNCE^) ^ ^^NCS)-^C2^^NCU) 
^ + i 5 ^q 

There is a barrier for NON adduct formation. The isotope ex¬ 


change rate is calculated to be very low ( [Gamallo et al.|2Q10| ) 
as the main exit channel is N 2 + O. Moreover the quintuplet 
surfaces are repulsive. We neglect this reaction. 


- + CO 


12 




This reaction was first mentioned by [Watson et al.|(|1976|) and 


has been experimentally studied in detail by 


Smith & Adams 


(1980) in the 80-500 K temperature range. Their data can be 
fitted ( Liszt|2007) ). We introduce a new formula allowing us 
to describe the full temperature range with a single formula. 

i^CO + HCO+ ^ CO + Ri^CO^ 

This react i on has also been studied experimentally by [Smith 


[& Adams] ( |1980[ ). The exothermicity of the reaction has been 


reported to be 9 K by [Langer et aL ( 1984| ) from theoretical 
studies whereas [Smith & Adams |T980|) proposed a value of 
12 + 5 K. However, a later study by Lohr ([1998[ ) leads to a 
value of 17.4 K. [Mladenovic & Roueff[ ( [2014[ ) have recon¬ 
sidered this reaction and have confirmed the value of ILohrl 
( 1998| ). We include this exothermicity value in the present 
work. Experimental points are fitted through a power law as 
given in Table 


- i^C^ + CN ^ i^C^ + i^CN 

There are no bimolecular exit channels for this reaction and 
two different pathways lead to isotopic exchange. In a di¬ 
rect reaction, NC ^ ^^CNC^ ^ ^^CN -f whereas 

the indirect pathway involves CN ^ ^^CCN^ ^ 

^^CN-f There is no barrier in the entrance valley for 
both cases at the M06-2X/cc-pVTZ level. Moreover, the 
cyclic transition state from CCN^ and CNC^ is located at 
-436kJ/mol below the reactant energy, leading to fast iso¬ 
merization. The corresponding capture rate constant is 3.82 
xio-9 X (r/300)-o-4o. 

- ^^C + CN ^ ^^C + ^^CN In a similar way to the ionic case, no 
exothermic bimolecular exit channels are available for this 
reaction and two reaction pathways occur for possible iso¬ 
topic exchange. NC ^ ^^CNC ^ ^^CN-f holds 
for the direct process and CN ^ ^^CCN ^ c-^^CNC 
^ ^^CN -f ^^C describe the indirect exchange mechanism. 
No barrier is found in the entrance channel and the transi¬ 
tion state from CCN to c-CCN is low (-298 
reactant energy), in good agreement with 
P002[ ) so that isomerization is expected to 
ture rate is computed as 3.0 x 10“^^ cm^ s“^ 

- + HCN ^ + H^^CN 

There are no exothermic bimolecular exit channels for this 
reaction and there is no barrier for HCNC formation but 
isotopic exchange requires isomerization through a TS lo¬ 
cated close to the reactant level, involving a TS located at 
-16 kJ/mol at the M06-2X/cc-pVQZ level but -f34 kJ/mol at 
the RCCSD(T)-F12-aug-cc-pVQZ level. Calculation of the 
rate constant for exchange is complex and similar to the ^^N 
-f CNC^ case. In the nominal model we neglect this reaction. 
We estimate nevertheless an upper limit for the rate constant 
by using a 1/T temperature dependence and assuming that 
the isotopic exchange rate constant balances the back disso¬ 
ciation at very low temperature (3K). 


kJ/mol below the 
Mebel & Kaiser 
be fast. The cap- 
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- 13c + C 2 ^ 

There are no exothermic product channels for this reaction 
and there is very likely to be no barrier. The capture rate 
constant is equal to 3.0 xl0“^^ cm^ s“^ 

- + CO ^ + CH 


This reaction could be an additional possibility to enhance 
as no barrier has been found in the entrance valley 
( Le Picard & Cano^|1998| ). The high pressure CH + CO 
^ HCCO association reaction rate constant is equal to 3 
xl0“^^ X (r/300)“^-^ between 53 and 294 K. The exchan ge 
rate has been measured to be ~ xl0“^^ ( Taatjes 1997| ) at 
room temperature. This value represents 1% of the associ¬ 
ation reaction rate constant at high pressure, which is ex¬ 
plained by a transition state localized at 6kJ/mole above the 
reactants energy as com puted at the M06-2X /cc-pVTZ level, 
in good agreement with |Sattelmeyer| ( [20041 ). We do not in¬ 
clude this reaction in our models given the high TS which 
should make this exchange process negligible at low temper¬ 
ature. 


Appendix B: ZPE values 

We revisit the ZPE values in light of several recent studies. 
We recall in Table the atomic masses of various isotopes 
which may be used to derive spectroscopic constants of isotopic 
molecules, as found in basic molecular spectroscopy textbooks 
( |Herzberg|1945|p^89] ). We only refer to the first order expansion 
terms for the purpose of computing ZPEs. 


Appendix B.0.1: Diatomic Molecules 

Eor diatomic molecules, energy levels are expressed as a Dun¬ 
ham expansion or equivalently as a sum of harmonic -l- an- 
harmonicity correction factors. The following expression is ob¬ 
tained for the ZPE, corresponding to v = 0. 

ZPE = -Tio + -T 20 = 2^0 - -^OJeXe (B.l) 

ojQ is the harmonic contribution of the vibrational frequency and 
ojeXe represents the anharmonic contribution. The reduced mass 
p dependence of ojq and ojeXe are 1/ and l/p respectively. 
The label computed in Table [^indicates the use of this property 
to compute the spectroscopic constants, in the absence of other 
information. 


Appendix B.0.2: Polyatomic molecules 

Eor polyatomic molecules, the following expression extends the 
diatomic formulae where the different vibrational degrees of 
freedom are included: 

Go = Z "'f Y (B-2) 

i i j 

(jL>i refers to the harmonic frequencies, di is the corresponding 
degeneracy and x/^ stands for the anharmonic terms. The sum is 
performed over the number of vibrational modes. In the case of 
polyatomic linear molecules, the number of vibrational modes 
is 3N-5 whereas it is 3N-6 in the general case where N is the 
number of nuclei in the molecule. 


Table 1. Atomic Masses in amu. from NIST 


Atom 

Mass 

H 

1.007825032 07 

D 

2.014 101 777 8 

N 

14.003074004 

15n 

15.000108 898 

C 

12.00 


13.003354 838 

0 

15.994 914 619 6 

180 

17.999 161 0 

e 

0.00054858 


Table 2. Spectroscopic constants of diatomic molecules in cm ^ and 
differences of ZPE with respect to the main isotopologue. 


Molecule 

N‘5n 
‘5N2 
NH 
‘^NH 
ND 
‘^ND 
NH+ 

‘5nh+ 

ND+ 

‘5nd+ 

CN 
‘^CN 
C'^N 
‘^NO 
‘^NO 
C 2 

‘^cc 

CO 

‘^co 

CH 
>3CH 
CD [14] 


2358.53 

14.30 

2319.01 

13.83 

2278.80 

13.35 

3282.2 

78.3 

3274.9 

77.95 

2399.0 

42.0 

2389.0 

41.6 

3047.58 

72.19 

3040.77 

71.87 

2226.93 

38.55 

2217.60 

38.23 

2068.68 

13.12 

2025.25 

12.57 

2036.70 

12.70 

1904.11 

14.08 

1870.1 

13.585 

1853.5 

12.755 

1817.3 

12.1 

2169.8 

13.3 

2121.4 

12.7 

2860.75 

64.44 

2852.15 

64.04 

2101.05 

34.7 


1175.7 

- 

1156.05 

28.3 

1136.06 

57.0 

1621.5 

- 

1617.9 

5.2 

1189.5 

621.6 

1183.6 

1505.7 

630 

1502.4 

4.8 

1103.8 

578 

1099.2 

585 

1031.1 

- 

1009.5 

31.1 

1015.2 

22.9 

948.5 

- 

931.7 

24.3 

923.98 

- 

905.61 

26.4 

1081.6 

- 

1057.5 

34.7 

1414.3 

- 

1410.1 

6.0 

1042.1 

535.5 


Ref oje ZPE A(ZPE) (K) 

[ 1 ] 

[ 1 ] 

[ 1 ] 

[ 2 ] 

computed 

[3] 

computed 

[4] 

[4] 

[4] 

computed 

[5] 

[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 

[ 10 ] 

[ 11 ] 

[ 11 ] 

[ 12 ] 

[13] 


[l]|Le Roy et al. 2006| [2] Ram & Bernath|2010 [3] Dore et al. 

2011 

[4] Colin|1989 [5] |Ram et al.||2010l [6] |Ram & Bernath 

2012 

[7] IColin & Bernath||20121 [8] | 

Henry et a 

L| 1978 [9] 

Danielak et al. 1997 [10] Zhang et al.| 

2011 [11] 

Guelachvili 

et al. 119831 [12] |Zachwieja||1995| [13] 

1 [Zachwieja 

119971 [14] 

Zachwieja et al.|2012 
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Table 3. Spectroscopic constants of triatomic molecules in cm ^ and 
differences of ZPE with respect to the main isotopologue. 


Molecule 

symmetry 

Ref 

(jJl 

1 : + 

(jJ2 

n 

(jJ2) 

ZPE 

A(ZPE) (K) 

HCN 

[1] 

3443.1 

727.0 

2127.4 

3512.25 

- 

H'^CN 

[1] 

3424.0 

720.6 

2093.0 

3478.6 

48.4 

HC'^N 

[1] 

3441.7 

725.9 

2094.0 

3493.75 

26.6 


[1] 

3422.6 

719.5 

2058.6 

3460.1 

75.0 

DCN 

[2] 

2702.5 

579.7 

1952.3 

2883.9 

904.1 

CCN 

[3] 

1967.2 

322.2 

1058.3 

1794.20 

- 

'^CCN 

[3] 




1778.19 

23.0 

C'^CN 

[3] 




1761.11 

47.6 

CC'^N 

[3] 




1775.65 

26.7 

HNC 

[4] 

3819.9 

463.8 

2064.3 

3369.9 

- 

HN'^C 

[4] 

- 

463.5 

- 



H'^NC 

[4] 

- 

461.5 

- 



N 2 H+ 

[5] 

3398.6 

695.5 

2297.7 

3508.6 

- 

'5NNH+ 

[5] 

3396.5 

694.4 

2259.3 

3487.5 

30.4 

N15nh+ 

[5] 

3382.4 

690.7 

2268.7 

3481.8 

38.5 

‘ 5 N 2 H+ 

[5] 

3380.7 

699.6 

2229.3 

3460.4 

69.3 

N 2 D+ 

[5] 

2719.8 

550.4 

2065.9 

2916.2 

852 

N‘5ND+ 

[5] 

2680.4 

544.4 

2058.8 

2888.4 

40.0 

'^NND-^ 

[5] 

2703.8 

549.0 

2041.0 

2895.2 

30.2 

CNC-" 

[6] 

1335 

165 

2040 

1852.5 

- 

C'5NC+ 

[6] 

1334 

161 

1997 

1826 

38.1 

HCO+ 

[7],[9] 

3229.6 

832.5 

2195.8 

3524.60 

- 

h‘3co+ 

[8],[9] 

3206.1 

824.9 

2160.6 

3488.24 

52.3 

DCO+ 

[7],[9] 

2645.9 

666.5 

1928.0 

2944.22 

835.1 


[l]|Maki et al.|(|2000l), [2] |Mollmann et al.|20Q21 [3] [Grant Hill| 


et al.|20111|Mitrushchenkov|2014i [4]|Maki & Mellau|2001| [5] 


Huang et al.|2010[ [61 Jensen & &aemer|1988| [7] [Martin et al.| 

1993 1 [8] private communication by T. Lee, based on second- 
order vibrational perturbation theory using a quartic force field 
determined at the CCSD(T)/cc-pVTZ level of electronic struc¬ 
ture theory, as in [7], [9] ZPEs are from |Mladenovic & Roueff| 

( [20T^ . 


Table 4. Spectroscopic constants of tetratomic nitrogen molecules in 
cm"^ and differences of ZPE with respect to the main isotopologue. 


Molecule Ref 

OJX 

(jJ2 

(j)2 

(jJ^ 

^^5 

ZPE A(ZPE) (K) 

symmetry 




n 

n 



HCNH+ 

[1] 

3484.2 

3187.3 

2159.9 

804.8 

649.4 

5987.3 

- 

H‘3CNH+ 

[1] 

3479.5 

3170.1 

2126.1 

799.8 

647.4 

5952.0 

50.8 

HC'^NH+ 

[2] 

3469.3 

3187.1 

2133.3 

804.7 

645.1 

5961.5 

37.1 

HCND+ 

[1] 

3221.9 

2681.1 

2025.9 

794.3 

513.1 

5364.2 

896.5 

DCNH+ 

[1] 

3474.0 

2601.0 

1934.8 

664.7 

622.8 

5388.0 

862.3 


Molecule 

Ref 


(jJ2 

(jJ2 

(jj/[ 

ZPE 

A(ZPE) (K) 

symmetry 


Ai 

Ai 

E 

E 



NH 3 

[3] 

3336.1 

932.4 

3443.6 

1626.3 

7429.2 

- 

'^NH3 

[3] 

3333.3 

928.5 

3435.1 

1623.2 

7413.0 

23.3 


[l]|Brites & Jutier|2012[ [2] Brites, private communication, [3] 
|Huang et al.|2Qll| 







































